Measuring Nonequilibrium Temperature of Forced Oscillators 
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The meaning of temperature in nonequilibrium thermodynamics is considered by using a forced 
harmonic oscillator in a heat bath, where we have two effective temperatures for the position and 
the momentum, respectively. We invent a concrete model of a thermometer to testify the validity 
of these different temperatures from the operational point of view. It is found that the measured 
temperature depends on a specific form of interaction between the system and a thermometer, which 
means the zeroth law of thermodynamics cannot be immediately extended to nonequilibrium cases. 



PACS numbers: 05.70.Ln, 05.40.Jc 



I. INTRODUCTION 



Temperature and entropy are basic concepts of thermodynamics which have clear definitions and meaning in equi- 
librium but which are not yet fully understood in nonequilibrium situations. In equilibrium thermodynamics, one way 
to introduce temperature is to define entropy somehow (e.g. through the adiabatic invariant) so that temperature 
can be introduced as a derivative of entropy with respect to energy; 



(i) 



O " dU' 

However, since nonequilibrium entropy has never been constructed in a consistent way, we cannot define nonequi- 
librium temperature in this manner. For example, we don't know whether the entropy is a measurable quantity in 
contrast to equilibrium cases where entropy difference between two states is measurable by heat produced in qua- 
sistatic processes. Furthermore, even if we can measure nonequilibrium entropy, we cannot obtain unique temperature 
£-»« | unless we properly set up the thermodynamic state space jl], ^, [| Q]; i.e. the value of the temperature depends on 
the choice of variables which we will fix through the differentiation of entropy with respect to energy jlj . 

Local equilibrium temperature, which we are familiar with, looses its validity for systems where the deviation from 
equilibrium ensemble is not negligible. Indeed, it is expected that equipartition of energy will no longer be valid, in 
such a way that different degrees of freedom may have different energy. For instance, some numerical simulations 
showed that nonequilibrium systems are anisotropic regarding with their kinetic energy || , which is never explained 

■ based on the local equilibrium assumption. 

So far, several authors have tried to seek the meaning of temperature beyond the local equilibrium picture utilizing 
I . the microscopic expression devised by Rugh || . 

> ■ 

where N is the number of degrees of freedom, m denotes mass of the microscopic particles, and <p is the interparticle 

■ potential. Note that we need arbitrary factors a and b to let the dimensions of the two terms (both of the numerator 
and of the denominator) be the same. We remark that, however, this arbitrariness has no influence on the value 
of temperature in equilibrium situations. Although Eq. (|^) is originally defined in microcanonical ensemble, Jepps 
et al. |7[ generalized this expression for canonical ensemble and presented it in a more general form. Furthermore, 
they applied it to numerical simulations of nonequilibrium stationary states in the presence of shear flow or heat flow. 
However, in nonequilibrium systems, we cannot have the unique value of temperature due to the arbitrariness of a and 
b || . This ambiguity seems quite natural since Rugh's expression is essentially based on equilibrium thermodynamic 
relation Eq. (|l|). Namely, the problem is carried over from the choice of variables to be fixed: we cannot reach the 
dynamical expression of nonequilibrium temperature unless the correct thermodynamic state space is set up. 

There is another way to define temperature which we call operational temperature in this paper. When finite closed 
systems are in contact, they finally equilibrate to have the same intensive quantity, which we identify with temperature 
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(i.e. the zeroth law of thermodynamics). Hence, it might be possible to measure nonequilibrium temperature by 
putting equilibrium thermometer in contact to nonequilibrium systems. A gedanken experiment has been proposed 
in order to clarify the meaning of nonequilibrium temperature from that point of view |J. In particular, Baranyai 
has performed numerical simulations on shear-flow or heat-flow systems in contact with thermometers and obtained 
some explicit values of operational temperature p]. However, since we don't have a theoretical framework in which 
the obtained values should be interpreted, those numerical data seem to be left alone. In other words, we cannot 
theoretically predict the value of operational temperature with which the numerical data should be compared. 

Since the meaning of temperature out of equilibrium seems to lack a sound theoretical basis, it is reasonable to pick 
up the simplest model in order to analyze theoretical aspects more easily. For that purpose, it seems that systems 
such as shear flow or heat flow are still excessively complicated. In this paper, we adopt a forced harmonic oscillator 
as a model system. Although it might be regarded as one of the simplest nonequilibrium systems, it is worth noting 
that time-averaged distribution functions of momentum and position are both Gaussian but characterized by different 
effective temperatures depending on the forcing frequency pOj . Thus, this system may provide a simpler example for 
the concept of temperatures than fluid systems under shear where different effective temperatures have been studied 
so far. We will simulate an experiment analogous to || , by letting a forced oscillator interact with another nonforced 
oscillator (in a different heat bath) which acts as a thermometer. Comparison of the respective results may be useful 
for clarification of the concept of temperature. 

The plan of the paper is as follows. First, some statistical properties of a forced harmonic oscillator in a thermal 
bath are recalled and interpreted in terms of nonequilibrium temperature. In the second section, in order to define 
temperature in a macroscopic point of view, a forced and an unforced oscillators situated in different thermal baths 
will be considered and the heat current between them will be calculated. In the final section we discuss the form of 
the entropy for some different choices of variables, and compare our result with those obtained by Baranyai in the 
framework of nonequilibrium molecular dynamics of fluid systems of soft spheres in shear flow. 



II. A MODEL SYSTEM: FORCED HARMONIC OSCILLATOR 



We assume that our model system is described by the following Langevin equation; 

x + ^x + n 2 x = Asmut + £(t) 7 (3) 

where the mass of the oscillator is taken as unit. The natural frequency of the harmonic oscillator is denoted by 
f2, and Asimot corresponds to external forcing. The noise term £(t) is assumed to be Gaussian white noise which 
satisfies 

(£(«)> = 0, (£(f)S(f))=2 7 /r 1 *(t-*'), (4) 

where (3 is the inverse temperature of the heat bath. 

Macroscopic or thermodynamic quantities should be defined by an appropriate averaging; usually ensemble-average 
or time-average. Throughout this paper, we will take time-averaged quantities as macroscopic variables, since the 
model system is periodic in time due to sinusoidal forcing. 

In our model, the external force gives power input to the oscillator, which may cause different influence on the 
average energy of momentum and of position. To have a physical idea of this influence, note that we have two kinds 
of relaxation times, each of which is related to position and momentum, respectively. We write the relaxation time 
of position as t x = 7f2~ 2 and the one of momentum as t p — 7 _1 . When the forcing period ta = 2ir/ui is longer 
enough than a relaxation time (i.e. ta 3> t x or r p ), the sinusoidal motion is averaged to yield a distribution function 
which deviates from the equilibrium one. In contrast, if the forcing period is comparable with (or shorter than) a 
relaxation time (i.e. ta < t x or t p ), the corresponding motion of position or momentum cannot follow the forcing and 
the distribution function is indistinguishable from the equilibrium one. For instance, if t p <^ta < t x , we may expect 
that the distribution function of position is not much changed from equilibrium, whereas the one of momentum is 
modified by the forcing. 

To see these circumstances explicitly, we will calculate the potential energy u con and the kinetic energy Ukin- Since 
the dynamics is linear, we can separate the ensemble-averaged motion and the fluctuation. After some textbook-like 
calculations we get 
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Ukin = ^7,+—^, (6) 
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where 



X Q = 



A 



^(fl 2 -LU 2 ) 2 + 7 2 w 2 ' 



The corresponding energy dissipation rate w is 



A 2 



w 



2 ( 7 w) 2 + (n 2 -u 2 ) 2 ' 

j(2u kin - (3- 1 ). 



(7) 

(8) 
(9) 



Then we wish to consider the distribution function described by the Kramers equation corresponding to our model 
Eq. (|). 



p(x,p;t) 
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p(x,p;t). 



The solution independent of the initial condition is 

P 



p(x,p; t) oc exp 



(p — ujX cos mty 



pn 2 



(x — Xq smutty 



(10) 



(11) 



where the time axis is shifted suitably. Since this solution corresponds to ensemble distribution at a given instant, in 
order to calculate time-averaged quantities, the distribution function Eq. (|ll| ) itself must be time-averaged, 



P(x,p) oc 



dtp(x,p;t), 



oc exp [- 1 (p 2 -nV) 



dt exp[— f3(pu>Xa coswt — CIxXq sinwt)]. 



(12) 



To perform this integral, we expand the integrand in Eq. (£3) up to second order in /3pujXo and /3Q 2 xXo (Gaussian 
approximation), providing that /3(wX ) 2 ^ 1 an d /3(^ 2 X )^5c 1. Namely, we assume that gain of the internal energy 
due to the external forcing is small compared with thermal energy. With this approximation, we get the following 
distribution function after simple calculations. 
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(13) 

(14) 
(15) 



We can see there are two kinds of temperature for one oscillator. Hereafter each temperature corresponding to position 
and momentum is called configurational temperature and kinetic temperature, respectively. 



III. OPERATIONAL TEMPERATURE 



The discussion in the previous section deals with the distribution function in the phase space and hence it may be 
rather microscopic consideration. From the thermodynamic point of view, the problem arises how the different mi- 
croscopic temperatures of Eqs. ([l4]) and (|l5| ) are connected to macroscopic measurements. Since those temperatures 
differ, the macroscopically measured temperature (i.e. operational temperature) may indicate different values depend- 
ing on the details of the connection between the system and the thermometer. In this section, we will investigate the 
problem by devising a concrete model for the temperature measurement. 

Here we will examine a situation which bears some similarities with the proposal by Jou and Casas- Vazquez ||. 
They defined a prototype of operational temperature in which two systems are in thermal contact: the one is kept 
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FIG. 1: Schematic of the numerical experiment for operational temperature. 



in a nonequilibrium steady state by means of heat flux, whereas the other is in equilibrium to act as a thermometer. 
In our setting, we consider two coupled oscillators in contact with different heat baths whose temperature can be 
controlled independently. The one is forced to stay away from equilibrium (the system) while the other is unforced to 
remain in equilibrium (the thermometer) . They are connected through a weak interaction potential V. The schematic 
picture of our situation is shown in Fig. |l} 

The dynamics of such oscillators will be written as 



x + jx + Sl 2 x + e 
a V + IV + e 



dV(x — y) 

O.i' 
dV(x ~ y) 

dy 



= Asm(ujt) 



(16) 
(17) 



where a denotes the mass of unforced oscillator. We let a and e to be small so that the disturbance of the forced 
oscillator by the thermometer will be weak. The noise terms £i(t) are again assumed to be Gaussian white noise but 
characterized by different temperatures denoted by /2T , i.e. 



(&(*)} = 0, (&(*)&(*')> = 2rf^%5(t - t'), 



(18) 



where is the Kronecker's delta (unity if i = j and zero otherwise). The similar system, in the absence of inertia 
and the external forcing, was studied in detail by Sekimoto Jn| . 

We will test two kinds of interaction terms; harmonic and bistable potentials. 
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The heat flux from the forced oscillator to the thermometer is evaluated as 

dV(x — y) . 

q = -e — 5 y- 

dy 



(19) 



(20) 



When both oscillators are left unforced (A = 0), the heat flux between both systems is proportional to the difference 
of the temperatures of the corresponding heat baths. When one of the oscillators is forced (A ^ 0), the unforced 
oscillator plays the role of thermometer. 

Our definition of the operational temperature is as follows. We fix the parameters of the forced oscillator (i.e. A, 
uj, (3\ and fi), and change the temperature of the heat bath for the thermometer (/?2)- There should be a certain 
value of @2 at which the average heat flux (q) vanishes. Then we identify the temperature of heat bath ft^ 1 with the 
operational temperature. Throughout the numerical simulations, we set 7 = 1.0, A = 1.0, j3\ = 1.0, e — 0.1, and 
a = 0.1. 

The results of the numerical simulations are shown in Figs. ^ and ^ for different parameters of the oscillator, 
where the values of (q) are rescaled suitably. Note that the zero-point of heat flux (equilibration point) is different 
depending on the interaction potential: thermometers indicate different values for the same system. Especially, the 
one with harmonic potential shows good agreement with the configurational temperature, while the other with bistable 
potential indicates a value close to the kinetic temperature. This tendency is unchanged for the system with different 
parameters at which the configurational temperature is higher than the kinetic one (Fig. ||), while the latter is higher 
than the former in Fig. |^. 
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FIG. 2: Average heat flux (q) as the function of the temperature of the heat bath of thermometer. The results of harmonic 
interaction are represented by x's, while the ones of bistable interaction are by +'s. The parameters are set as ui = 0.25 
and £1 = 1.0, where the corresponding kinetic and configurational temperatures are 1.08 and 1.49, respectively (shown by the 
arrows) . 
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FIG. 3: The same graph with Fig. g but with different parameters such that uo = 1.0 and Q. — 0.1, where the corresponding 
kinetic and configurational temperatures are 1.37 and 1.04, respectively (shown by the arrows). 



In short, the operational temperature of a thermometer with harmonic potential is almost the configurational 
temperature, and the one with bistable potential is close to the kinetic temperature. It is concluded that various 
interactions show various temperatures whose values are ranged between kinetic and configurational temperatures. 
Therefore, it is plausible that the zeroth law of thermodynamics should be formulated only in a very restrictive form 
for nonequilibrium cases. 



IV. CONCLUDING REMARKS 



In this concluding section, we will pay some attention to entropy corresponding to the nonequilibrium system, and 
to a comparison with the results obtained by Baranyai and coworkers by means of the nonequilibrium molecular 
dynamics. 



A. Entropy for nonequilibrium steady states 

In this subsection we will see some possible extended Gibbs relations of forced harmonic oscillators. The motivation 
is to discuss temperature in thermodynamic point of view and to compare it with the operational temperature obtained 
in the previous section. The distribution function Eq. ([l^) based on the Gaussian approximation may be considered 
to give the maximum entropy with the constraint that the second moments of momentum and position are given. 
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The corresponding Gibbs relation will thus be 

ds = (i con du con + /3 k inrfukin, (21) 

where s denotes entropy per oscillator. This may be a natural extension of the equilibrium Gibbs relation by consid- 
ering that each degree of freedom is a thermodynamic system by itself with different temperature. This simple form 
of the Gibbs relation is due to the Gaussian approximation of the distribution function. If they are not Gaussian, 
we would need additional independent variables; e.g. higher order moments. Recall, anyway, that systems with two 
temperature are common in noncquilibrium physics, for example in plasmas (where electrons and ions may have 
different temperatures), or in metals or semiconductors (where electrons may exhibit a temperature different from 
that of the lattice) . 

However, Eq. (|TJ) is not the only candidate. We can consider arbitrary linear transformation of it con and Ukin such 
that 

11 = Ucon + Ukin, (22) 

y = aiu con + a 2 u kin , (23) 
where a\ ^ a 2 - By using these new variables, we can rewrite Eq. (|l]) as 

^l/^kin /^con /^kin (24) 

Since u is the total energy of the oscillator, it may be possible to define temperature analogous to Eq. (|l]); 
which yields 

q _ Q2/3con - OlAdn 

a 2 — ot\ 



(26) 



Here we use another notation 9 as the (inverse) temperature of the nonequilibrium system, which clearly depends on 
the choice of ot\ and ot 2 . Namely, the thermodynamic temperature 9 depends on the choice of the new variable y. 

Although the nonuniqueness of temperature has been argued, e.g. in the context of extended irreversible thermo- 
dynamics |Q] , the criteria for the choice of new variables is still unknown. At least there are some necessary conditions 
for the choice of y: 



• it must be extensive, and macroscopically observable. 

• the entropy must be convex regarding with the new extensive variable y. 



In Eq. (21), the convexity is identical with that f3 con and /3kin are nonincreasing functions of u con and itkin , respectively. 
This is obvious since /3 con = (2u con ) -1 and /3ki n = (2ukin) _1 - Of course there are other choices satisfying convexity. 
Say we set a± = 1 and a 2 — 0, then 

ds = Pkindu + (/3 con - /3 k in)rfWcon, (27) 

where the thermodynamic temperature coincides with the kinetic temperature. 
Another possibility is to adopt the entropy production rate a 

a = /3u; = 7(2/3ukin-l), (28) 
where Eq. @ is recalled. Since da — 2"fj3du^\ n , from Eq. ( ^l| ) we get 

ds = (3 C ondu - ^ con ^ km da : (29) 

where the thermodynamic temperature becomes the configurational temperature. Using Eqs. (|l4|) and (|l5|), we can 



further rewrite Eq. (29) as 

2 7 2 V' w 



ds = f3 con du - ( 1 - — ) ada, (30) 
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where we can see the second-order contribution of the flux a to the entropy. Note that the above expression with 
energy and entropy production rate are analogous to that in extended irreversible thermodynamics |Lj, where the 
usual thermodynamic variables and the fluxes are taken as independent variables. In this case we have /3 con as the 
thermodynamic temperature. In addition, we remark that the second terms of Eqs. (p7|) and (|29| ) vanish at resonance 
where w = Q. It is identical with 



dy 

Also the entropy production rate a is maximum at resonance. 

As we have seen so far, thermodynamic temperature defined through extended Gibbs relation depends on the 
choice of the new variable. We have expected that the operational temperature would be the criterion for choosing 
the new variable, which was one of the motivations of our study. However, in the previous section, we have seen that 
different thermometers read different temperatures, which means impossibility of defining the unique nonequilibrium 
temperature even in this simple model system consisting of one degrees of freedom. The absence of a unique operational 
temperature can be a serious problem for the construction of thermodynamics: at least the formulation of the zeroth 
law is not immediate and, if possible, would be a very restrictive form in contrast to equilibrium cases. Sasa and 
Tasaki have already pointed out this kind of operational restriction which results from the anisotropy of pressure in 
a macroscopic heat conducting system |12j . 



B. Comparison with the results of nonequilibrium molecular dynamics 

In this paper we have examined kinetic, configurational, and operational temperatures in a forced harmonic oscilla- 
tor. As was mentioned in the introduction, similar situations have been examined by Baranyai |^| based on the same 
motivation. While he studied systems consisting of soft spheres under shear flow or in the presence of heat current 
using techniques of nonequilibrium molecular dynamics, the system analyzed here is much simpler than that. It must 
be noted that Baranyai used the Nose-Hoover type dynamics which removes energy from the system as dissipation, 
whereas we adopt the Langevin equation to represent the effect of heat baths. Despite these differences, we think 
that it is still worth comparing these results since we are looking for general thermodynamic concepts which should 
be largely independent of the microscopic details of the system. 

The definitions of temperature in the works of Baranyai is based on the Rugh's microscopic expression of Eq. (|^). 
As we have mentioned in the introduction, the expression itself cannot define temperature uniquely in nonequilibrium 
states due to arbitrary factors a and b. Instead of the unique temperature, again we have two kinds of temperature; 
configurational and kinetic temperatures which can be defined without the arbitrariness. 

Pcon = {ww)' (32) 

N 

Adn = . m ^y (33) 
(Si m / 

Theses expressions yield /? con = (2u con ) _1 and /3km = (2itkin) _1 which coincide with our results obtained by Gaussian 
approximation. Baranyai calculated these expression for temperatures in the mentioned fluid system of soft spheres at 
various values of the shear rate, and found that the configurational temperature is higher than the kinetic temperature, 
whereas the situation is opposite in systems with charge current. Furthermore, these temperatures turned out to be 
anisotropic: i.e. they take different values for different spatial directions. In our situation, the relation between the 
configurational and the kinetic temperatures depends on the ratio Q/u). In general, as Baranyai has discussed, they 
will depend on the characteristics of the system and the external forcing responsible for the nonequilibrium situation. 

In addition, Baranyai has proposed an operational temperature by devising the concrete model which emulates a 
physical thermometer in contact with the fluid. The thermometer consists of the same particles as the fluid's, but do 
not feel the effect of shear flow nor the thermostatting: they interact only with the fluid particles. This thermometer 
seems to read definite values of temperature regardless of the mass and the number of thermometer particles. However, 
the dependency of interaction potential between the fluid and the thermometer is not discussed. Taking our result 
into consideration, the operational temperature will depend on the interaction between the system and thermometers. 
Indeed, Hoover et al. have discussed the ideal gas thermometer which reads the kinetic temperature while 
Baranyai's thermometer reads the value which is closer to the configurational temperature. However, since Hoover 
et al. ignored the anisotropy of the kinetic energy, it is not apparent what value the ideal gas thermometer reads 
when it is actually in contact with a nonequilibrium system. More numerical simulations and real experiments 
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on the operational temperature are needed for the clarification of nonequilibrium temperature and the underlying 
nonequilibrium thermodynamics. 

In summary, we have found that the operational temperature depends on the details of interaction between the 
system and a thermometer, which is never seen in equilibrium situations. The fact may shove a strong restriction on 
the extension of the zeroth law for nonequilibrium systems. Of course, as our analysis is confined to one-dimensional 
systems, there may arise another problem for two- or three- dimensional systems regarding with the relation between 
the anisotropy of temperatures and the operational temperature. In addition, the Gaussian approximation of the 
distribution functions make especially easy to define temperatures, while additional conceptual problems would appear 
if the distribution function deviates far from the Gaussian. For instance, there have been some maximum-entropy 
analyses of nonequilibrium radiation, where nonequilibrium temperature or quasitemperature per mode have been 
defined in terms of the nonequilibrium populations of the different modes, in the context of a generalized Planck 
statistics instead of a classical Boltzmann statistics The present situation has the advantages of the higher 

simplicity and of the possibility of devising numerical simulations concerning operational temperature which have not 
been done before. 
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